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A new version of the invariant imbedding theory for the propagation of coupled waves in inho- 
mogeneous media is applied to the mode conversion of high frequency electromagnetic waves into 
electrostatic modes in cold, magnetized and stratified plasmas. The cases where the external mag- 
netic field is applied perpendicularly to the direction of inhomogeneity and the electron density 
profile is linear are considered. Extensive and numerically exact results for the mode conversion 
coefficients, the reflectances and the wave electric and magnetic field profiles inside the inhomoge- 
neous plasma are obtained. The dependences of mode conversion phenomena on the magnitude of 
the external magnetic field, the incident angle and the wave frequency are explored in detail. 



I. INTRODUCTION 

Mode conversion from one type of wave mode to an- 
other type of mode at resonance points in an inhomoge- 
neous plasma is a very important phenomenon in various 
branches of plasma physics^2*2i4 In a previous paper^ 
we embarked on a project of developing an exact theory 
of mode conversion in inhomogeneous plasmas, which is 
based on a new version of the invariant imbedding the- 
ory developed by us earlier In that paper, we consid- 
ered the mode conversion of the simplest kind, that is 
the conversion of obliquely-incident p-polarized electro- 
magnetic waves into electrostatic plasma oscillations in 
cold, unmagnetized and stratified plasmas pSuLifiiil Using 
the invariant imbedding method ^iSiiiiiiSiiSiii we were 
able to calculate the mode conversion coefficient and the 
electromagnetic field distribution for linear and parabolic 
plasma density profiles in a numerically exact manner. In 
some cases, there were substantial discrepancies between 
our results and previous theoretical calculations. 

The mode conversion phenomena in cold, unmagne- 
tized plasmas are described by a single wave equation. 
A large number of wave propagation problems in more 
general situations, however, are described by two or more 
coupled wave equations. For the study of those problems, 
the invariant imbedding method needs to be generalized. 
This generalization was already made in Ref. 7 and has 
been successfully applied to various problems where sev- 
eral coupled waves propagate in inhomogeneous media. 
The propagation of electromagnetic waves in magnetized 
plasmas also involves several coupled waves. In this pa- 
per, we apply the theory developed in Ref. 7 to the mode 
conversion of high frequency electromagnetic waves into 
electrostatic modes in cold, magnetized and stratified 
plasmas. 

Even in the relatively simple cases where the plasma 
is cold and stratified, the mode conversion of electro- 
magnetic waves in magnetized plasmas is a highly com- 
plex problem. There are three independent directions in 
three-dimensional space, which are the directions of inho- 



mogeneity, the external magnetic field and the incident 
wave vector. The wave propagation and mode conversion 
phenomena depend on the relative orientations of these 
three directions. In the present work, we limit our at- 
tention to one special configuration where the external 
magnetic field is applied perpendicularly to the direc- 
tion of inhomogeneity and plane waves propagate in the 
plane spanned by the magnetic field and the inhomogene- 
ity direction. Other configurations will be considered in 
forthcoming papers. 

The theory developed in this paper can be applied eas- 
ily to multi-species plasmas, which consist of electrons 
and different kinds of heavy ions. In this paper, we will 
consider only the propagation of high frequency waves 
and assume heavy ions to be immobile. Rather compli- 
cated, but very interesting mode conversion phenomena 
associated with low frequency waves and ion density os- 
cillations will be subjects of our future publications. 

Unlike in the unmagnetized case, the wave propaga- 
tion in magnetized plasmas is anisotropic. In a uniform 
plasma where the wave vector is not precisely parallel to 
the external magnetic field, the ordinary (O) wave and 
the extraordinary (X) wave are two independent eigen- 
modes. When the wave frequency u> is bigger than the 
electron cyclotron frequency uj c , there can appear the O 
and X mode cutoffs and the X mode resonance inside 
the plasma. On the other hand, if the wave frequency 
is smaller than oj c , the O mode resonance can appear as 
well as the O and X mode cutoffs. In inhomogeneous 
plasmas, the O and X modes are no longer true eigen- 
modes and are coupled to each other. These modes can 
also be converted to upper hybrid oscillations at the res- 
onances. The presence of both a strong coupling between 
different modes and the singularity associated with reso- 
nance points makes the solution of the mode conversion 
problem very difficult. As will be demonstrated later, our 
theory allows us to obtain an exact solution of this diffi- 
cult problem in a remarkably simple and efficient manner. 

Mode conversion in magnetized plasmas 
has been studied extensively over several 
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We will not attempt to give a detailed review of the 
literature here. Many of the previous works have 
used one approximation scheme or another, most 
frequently the WKB method. In a series of papers, 
Mj0lhus calculated the mode conversion coefficient for 
cold, magnetized plasmas primarily using the WKB 
method^SL 2 ^ This work, which was reviewed in 
Ref. 4 together with theoretical results of many other 
researchers, often predicts a complete mode conversion 
for certain parameters, in great contrast with our 
results. This discrepancy might be explained in large 
part by one crucial difference between Mj0lhus' and our 
calculations. In the present paper, we consider only the 
cases where linearly-polarized waves are incident from 
a vacuum onto a stratified plasma, whereas Mj0lhus 
and other researchers have considered the cases where 
O or X waves are incident from a uniform plasma. 
Upon entering the plasma, linearly-polarized waves will 
be split into the O and X mode components. In the 
case of perpendicular inhomogeneity considered in this 
paper and when the waves are incident on a stratified 
plasma normally, mode conversion is caused solely by 
the X mode component. In this case, we find a perfect 
agreement between our result for incident s waves shown 
in Fig. ^and Mj0lhus' result for incident X waves. When 
the waves are incident obliquely and for sufficiently large 
values of the parameter w defined in Sec. IIV Al mode 
conversion is dominated by the incident O mode. This 
case is closely related to the theory of radio windows in 
the ionosphere and magnetosphere^^* 4 * 2 ^ 3 . We find that 
the critical parameter values where the mode conversion 
coefficient takes the maximum value agrees extremely 
well with the analytical formula, Eq. (|28[) . obtained 
by previous theories, for sufficiently large w. This is 
another evidence for the consistency between our theory 
and previous theories. For quantitative comparisons 
with previous results, we need to generalize our method 
to the cases where O or X waves are incident from a 
uniform plasma. This generalization will be pursued in 
future research. Since our theory is capable of producing 
exact results for realistic plasma density profiles, we 
expect it will be useful in the interpretation of various 
radio modification experiments. 

In Sec. |nl we present the details of a new version of 
the invariant imbedding theory for the propagation of 
coupled waves in stratified media. In Sec. IIII Al we in- 
troduce the coupled wave equations for electromagnetic 
waves in cold, magnetized plasmas with perpendicular in- 
homogeneity. In Sec. IIII 51 we describe two linear plasma 
density profiles used for the calculations in this paper. 
The results of our calculations are presented in Sec. IIVI 
In Sec. IIV Al we study the magnetic field dependence 
of mode conversion coefficients. 



In Sec. IIVBI we con- 
sider the dependence of mode conversion coefficients on 
the incident angle in detail. We also consider the spatial 
dependence of electric and magnetic field intensities. In 
Sec. IIV CI we discuss the frequency dependence of mode 



conversion coefficients. We conclude the paper in Sec.lVl 
In Appendix^] we give a detailed derivation of Eq. l|19f) . 



II. INVARIANT IMBEDDING THEORY OF 
THE PROPAGATION OF COUPLED WAVES IN 
STRATIFIED MEDIA 

In Ref. 7, we have developed a new version of the 
invariant imbedding theory for the propagation of cou- 
pled waves in arbitrarily-inhomogeneous stratified media. 
Since this theory will play a central role in our develop- 
ment of the theory of mode conversion, we give a detailed 
presentation of it here. 

We consider N coupled waves propagating in a strat- 
ified medium, where all physical parameters depend on 
only one spatial coordinate. We choose this coordinate 
as the z axis and assume the inhomogencous medium of 
thickness L lies in < z < L. We also assume that all 
waves propagate in the xz plane. The x component of 
the wave vector, q, is then a constant. In a large class 
of interesting problems, the wave equation of N coupled 
waves in stratified media has the form 

^ " f^W^ + [£(z)K 2 M(z) - q 2 l] ^ = 0, (1) 



dz 



dz 



where ip = (ipi , • • • , iPn) T is a vector wave function and £ 
and M. are N x N matrix functions. We assume that the 
waves are incident from a vacuum region where z > L 
and transmitted to another vacuum region where z < 0. 
I is a unit matrix and K is a diagonal matrix such that 
Kij = kiSij, where ki is the magnitude of the vacuum 
wave vector for the i-th wave. In the problem considered 
in this paper, all fcj's are equal. However, there exist 
many cases where fcj's take different values. It is useful 
to rewrite Eq. (JIJ as 



d 

dz 



£-\z) 



dip 
dz 



+ [K 2 M(z)-q 2 £- 1 {z)]ip = 0. (2) 



From this form, we can see that £~ l dip / dz is continuous 
everywhere except at some singular points. 

Eq. is suitable when the matrices £ and M. do 
not contain the constant q. In the cases where the q 
dependence cannot be separated as simply as in Eq. QJ, 
it is more convenient to write the wave equation as 



d 2 ip _ 

dz 2 dz dz 



£(z)K 2 M{z)ip = 



(3) 



and define £ and M. as g-dependent quantities. By as- 
signing £(z) and M.(z) suitably, we are able to describe 
various kinds of waves in a number of stratified media us- 
ing either Eq. fl} or 0. A wide variety of mode conver- 
sion phenomena observed in space and laboratory plas- 
mas can be studied using these equations. 

We generalize Eqs. and J3J) slightly, by replacing the 
vector wave function tp by an N x N matrix wave func- 
tion 'J', the j-th column vector ■ ■ ■ ,^Nj) T of which 
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represents the wave function when the incident wave con- 
sists only of the j-th wave of unit amplitude. Our main 
interest is in the reflection and transmission coefficient 
matrices r — r(L) and t = t(L). The wave function in 
the incident and transmissive regions are expressed in 
terms of r and t: 



iP(L-z) + e iP(z-L) ^ z > L 



e- lPz t, 



z < ' 



(4) 



where P is a diagonal matrix satisfying Py = piSij and pi 
is the negative z component of the vacuum wave vector 
for the i-th wave. It is important to write the products of 
two matrices, e lP ( z-i ) r and e~ lPz t, in the given order. 
From Eq. QJ, we obtain = I + r and *(0) = t. 

Let us introduce a matrix 



T exp [i f* dz" £{z")P] , z> z' 
t exp [-//* dz" £{z")P] , z < z' 



(5) 



where T and T are the time-ordering and anti-time- 



ordering operators respectively. When applied to a prod- 
uct of matrices which are functions of z, T (T) ar- 
ranges the matrices in the order of decreasing (increas- 
ing) z. For instance, T [£ {z\)£ (z 2 )] = £{z%)£{z\) and 
f [£{z 2 )£{zi)\ = £{z 1 )£(z 2 ), if z 2 > zx. Though, strictly 
speaking, T has to be called a space-ordering operator, 
we keep using the standard terminology, time-ordering 
operator, commonly used in quantum mechanics^ Ex- 
panding the matrix exponentials in Eq. in a Taylor 
series and using the definitions of T and T, we can easily 
prove that g(z, z 1 ) satisfies the differential equations 



0_ 

dz 
d_ 

dz 



g(z,z') = i sgn(z - z') £(z)Pg(z, z'), 
7 g(z,z') = -isga(z-z') g{z,z')£(z')P. (6) 



Using Eqs. JSJ and 10, the wave equation Q is trans- 
formed exactly to an integral equation 



9(z,L)=g(z,L) 



dz'g{ Zl z') [£{z')P - PM(z') - q 2 p- 1 M{z') + q 2 P~ 1 £- 1 (z')] L), 



r 



(7) 



where we have considered $ as a function of both z and where 
L. We replace L by a variable / and take a derivative of 
this equation with respect to I. Then we obtain 



&iF(z,l) 
dl 



i^(z,l)a(l) + ®(z,l), 



(8) 



t(l) = £{l)P - - [£{l)P - PM{1) - g 2 r'M(i) + q 2 p- 1 £- 1 {l)] I), 



dz'g{z, z') [£(z')P - PM{z') - q 2 p- 1 M{z l ) + q 2 P~ l £~ l {z')\ $(z', I) 



(9) 



We notice that $(z, I) satisfies an equation similar to 
Eq. |0 except that there is no source term (that is, 
g(z, I)). This implies Q(z, I) = and then we have 



&S>(z,l) 
dl 



M(z,l)a(l). 



(10) 



Taking now the derivative of ^(l, I) with respect to I, we 
obtain 



dV(l,l) _ &9(z,l) 



dl 



dz 



dV(z,l) 



z=l 



dl 



i£{l)P[r{l)- I]+W{l,l)a(l), (11) 



where we have used Eq. (|10[1 and the fact that 
£~ 1 {z) d^(z, l)/dz is continuous at z = I. Since ^(l, I) — 
I + r(l), we find the invariant imbedding equation satis- 
fied by r(l): 
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J = i[r(l)£(l)P + £(l)Pr(l)] 

~\W) + I] [SQ)P - PM(l) - q 2 P- 1 M{l) + q 2 p- l £- 1 {l)] [r(0 + /]. (12) 



I 

Similarly by setting z = in Eq. (JTUJ, wc find the invari- ant imbedding equation for t(l) (= ^(0, 1)): 
I 



-~t(0 [£ (OP - PM(l) - q 2 P~ 1 M(l) + (fP^e- 1 ^)] [r(J) + /]. (13) 



These equations are supplemented with the initial condi- 
tions, r(0) = and t(Q) = I. 

We integrate the coupled differential equations i|12|) 
and (| 1 3|) numerically from I = to I — L using the ini- 
tial conditions and obtain the reflection and transmission 



coefficient matrices r and t as functions of L. The invari- 
ant imbedding method can also be used in calculating the 
field amplitude ^(z) inside the inhomogeneous medium. 
Rewriting Eq. (TDJ, we get 



dV(z,l) 
dl 



M(z,l)£(l)P 

I) [£{l)P - PM(l) - q 2 p- 1 M(l) + q 2 p- 1 £- 1 (l)] [r(l) + I}. 



(14) 



For a given z (0 < z < L), the field amplitude ^(z, L) is 
obtained by integrating this equation from I — z to / = L 
using the initial condition ^{z, z) = I + r(z). 

Our invariant imbedding equations, Eqs. (|12I14|) . have 
been applied successfully to the study of the propagation 
of electromagnetic waves in multilayer systems made of 
chiral media and of phase-conjugating media. 7 



III. COUPLED WAVE EQUATIONS FOR COLD, 
MAGNETIZED PLASMAS WITH 
PERPENDICULAR INHOMOGENEITY 

A. Coupled wave equations 

Eqs. I|12|l . (|13|l and (|14|l are the starting point in our 
exact analysis of a variety of wave coupling and mode 
conversion phenomena. In the rest of this paper, we ap- 
ply these equations to the high frequency wave propaga- 
tion and mode conversion in cold, magnetized plasmas. 
For the sake of simplicity, we restrict our attention to the 
cases where the external magnetic field is applied perpen- 



dicularly to the direction of inhomogeneity. 

We assume that the plasma density varies only in 
the z direction and the uniform external magnetic field 
Bo (= BqSc) is directed in the x direction. The cold 
plasma dielectric tensor, e, for high frequency waves in 
the present geometry is written as 



(15) 











ic-2 


-«£2 


(~1 



where 



C2 



£3 



LO 2 (w + if) 



(oj + iv) — u>1 



1 



(lu + iv) — iJ], 

Up 

uj (u) + iv) 



(16) 



We note that the indices 1, 2 and 3 in e 1; e 2 and €3 have 
no direct relationship with the spatial coordinates x, y 
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and z. The constant v is the phenomenological collision 
frequency. In Gaussian units, the electron plasma fre- 
quency Up and the electron cyclotron frequency uj c are 
given by 



47re 2 . . eB 

n(z), oj c = , 

m P m„c 



(17) 



where m e and — e are the electron mass and charge re- 
spectively. The spatial inhomogeneity of plasmas enters 
through the z dependence of the electron number density 
n. 

For monochromatic waves of frequency to, the wave 
equations satisfied by the electric and magnetic fields in 
cold, magnetized plasmas have the form 



-V x (V x E) + -5-e • E = 0, 



-V x (e _1 • V x B) 



-B = 0. 



(18) 



In this paper, we restrict our interest to the cases where 
plane waves propagate parallel to the xz plane, which is 
sometimes called the magnetic meridian plane. We as- 
sume that all field components depend on x and t through 
a factor e l [( kBinS ) x -ut\ ^ w ]j ere ^ = ^jc and 9 is the inci- 
dent angle, and have no y dependence. In this situation, 
we can eliminate E Xl E z , B x and B z from Eq. I|18|) and 
obtain two coupled wave equations satisfied by the z- 
dependent field amplitudes E y — E y (z) and B y = B y (z), 
which turn out to have precisely the same form as Eq. (jSJ) 
with 



K=-I, 

c 







)■ 


£3 




£l - 


- 2 

ZL- - s in 


ei 




i^sinfl 







1 



sin 2 



(19) 



A derivation of this crucial result is given in Appendix 
A. 

We use Eqs. {T21, {T3J) and in calculating the re- 
flection and transmission coefficients. In our notation, 
Tn (^21) is the reflection coefficient when the incident 
wave is E y (that is, s-polarizcd) and the reflected wave 
is E y (By). Similarly, r 2 2 (^12) is the reflection coeffi- 
cient when the incident wave is B y (that is, p-polarized) 
and the reflected wave is B y (E y ). Similar definitions 
are applied to the transmission coefficients. The re- 
flectances and transmittances are defined by R L j — \rij\ 2 
and Tij = \Uj\ 2 . When the dielectric permittivities of 
the incident region and the transmissive region are the 
same, we can calculate the wave absorption by A, = 
1 — R\j — i?2j — T\j — T2j (j = 1, 2). If a mode conversion 
occurs, this quantity is nonzero even in the limit where 
the damping constant v goes to zero£ In that limit, we 
will call Aj as the mode conversion coefficient. 



Finally, we use Eqs. (|14f> and l|19(l in calculating the 
electromagnetic field amplitudes E y (z) and B y (z) inside 
the inhomogeneous region. Once we obtain E y and B y , 
we can use Eq. (IA2|) to calculate other field components 
E x , E z , B x and B z . 



B. Plasma density profiles 

Our theory can be applied to any function n(z) for the 
electron density profile. In the present work, we confine 
our calculations to two kinds of linear profiles. In the 
first one, we assume that the density is given by 



if z > L 

n(z) = \ n ^- if < z < L 
if z < 



(20) 



where the constant A is the scale length for the linear 
density profile, and the wave frequency is fixed at the 
value of the local plasma frequency at z = L — A, that is 
u) = (47rnoe 2 /m e ) 1 / 2 . In this case, the dielectric tensor 
for < z < L is given by Eq. 115|) with 



L-z 



1 



(■2 



£3 



A 



(1 



■ ~\2 
W) 



L-z 



Y 



1 



A (H 

L-z 



ivf - Y 2 ' 



A 1 



(21) 



where v = v/lo and Y = lu c /uj. For z > L or z < 0, 
the dielectric tensor is equal to a unit matrix. We will 
use this tensor in studying the dependencies of various 
quantities on the incident angle and the strength of the 
external magnetic field. 

We are also interested in the frequency dependence 
of mode conversion phenomena. For this purpose, we 
consider a linear density profile for < z < L such that 
the plasma frequency satisfies 

1/2 



■j p = 2nf 



1 + 24 



1 

L 



(22) 



where fo is the local (nonangular) plasma frequency at 
z = L measured in MHz. We introduce nonangular fre- 
quencies 

f =2^' /c =27' U =2^' (23) 

also measured in MHz. Then, for < z < L, the dielec- 
tric tensor is given by Eq. I|15(l with 

fo (/ + ifu) 



£1 = 1- 



/ 



C2 



£3 



fo fc 



1 + 24 



('-1) 



(/ + KUf-.fi 2 
Po 

f{f + iU) 



1 + 24 



1-- 

L 



1 + 24 



z 

I - - 



(24) 
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If z > L or z < 0, the dielectric tensor is equal to a unit 
matrix. For specific calculations, we will assume fa = 1 
MHz and f c — 1.5 MHz. Then the local plasma frequency 
varies from 1 MHz at z = L to 5 MHz at z = 0. 



IV. RESULTS 

A. Magnetic field dependence 

In this subsection, we study the dependence of mode 
conversion coefficients on the magnitude of the external 
magnetic field Bq, or equivalently on the cyclotron fre- 
quency lo c . We first consider the simplest case where 
linearly-polarized plane waves are incident normally (i.e. 
9 = 0) on a layer of inhomogcncous plasma with the lin- 
ear density profile given by Eq. (|20[1 . In this case, it is 
well-known that there is a resonance associated with X 
waves if the wave frequency is equal to the upper hybrid 
frequency uj u h = (t^ + w 2 ) 1 ^ 2 , whereas there is no reso- 
nance associated with O waves. For L > A, the resonance 
condition u> = uj u h can be satisfied if Y < 1 (i.e. lj > oj c ). 
The coordinate of the resonance point zr is given by 



z R = L - A (1 - Y 2 



(25) 



Since we are considering the cases where waves are inci- 
dent from a vacuum region, we assume that the incident 
wave is either s- or p-polarized. In the case of normal 
incidence, p waves are the same as O waves and s waves 
are similar to X waves. Therefore we expect that there is 
no mode conversion when p waves are incident normally, 
whereas some portion of normally-incident s waves will 
convert to upper hybrid oscillations. 

Actual calculations using our invariant imbedding the- 
ory show that the p wave absorption A2 vanishes indeed 
when the damping constant v is very small. In Fig. ^ 
we plot the s wave mode conversion coefficient A\ as a 
function of the dimensionless parameter w = C t 1 ^\/Y 
(oc ^/Bo) f° r several values of C = wA/c and for 8 = 0, 
v = 10~ 8 and L/A = 20. We notice that A\ approaches 
a universal curve as £ increases to large values. We find 
this universal curve agrees remarkably well with Fig. 4 
presented in Ref. 4. 

Next we consider the cases where linearly-polarized 
waves are incident obliquely at an angle 9. In those cases, 
both s and p waves can generate X mode components in- 
side the inhomogeneous plasma and be converted to up- 
per hybrid oscillations at z — zr. In Fig. [3 we show the 
s (p) wave mode conversion coefficient A\ (A2) as a func- 
tion of w for several values of the parameter v = C 1 / 3 sin# 
and for ( = 500, v = 10~ 8 and L/A = 20. As the incident 
angle increases, A2 grows from zero to finite values. 

In Fig.|SJ we plot the reflectances i?n, R22 and R12 for 
C = 500, v = 10~ 8 and L/A = 20. The parameter v is 
equal to 1 in Fig. 3(a) and 2 in Fig. 3(b). It can be proved 
that R21 is always identical to R\2- We observe that in 
the parameter region where there is no mode conversion, 
the reflectances are rapidly oscillating functions of w. 



B. Incident angle dependence and the 
electromagnetic field distribution 

In this subsection, we consider the dependence of mode 
conversion coefficients on the incident angle 9 in detail. 
We also consider the spatial dependence of electric and 
magnetic field intensities associated with the electromag- 
netic wave. In Fig. 0] we show the s (p) wave mode con- 
version coefficient A\ (A2) as a function of the parameter 
v, which is proportional to sin0, for several values of w 
and for ( = 500, v = 10~ 8 and L/A = 20. In the absence 
of the external magnetic field, the s wave mode conver- 
sion coefficient A\ vanishes, whereas the p wave mode 
conversion coefficient A2 agrees with the result for un- 
magnetized cases reported in Ref. 7. When Bq > 0, A\ 
becomes nonzero. For w values well over 1, mode con- 
version occurs in a narrow range of v or 9. This range 
is called the radio window and becomes narrower as w 
or £ increases. In Fig. 5, we plot the reflectances Rn, 
R22 and R12 for £ = 500, v = 10~ 8 and L/A = 20. The 
parameter w is equal to 0.5 in Fig. 5(a), 1 in Fig. 5(b) 
and 2 in Fig. 5(c). 

The occurrence of radio windows is closely related to 
the coupling between the O mode and the X mode near 
the position where u> — uj p , which has the coordinate 
z = L — A (= zq) in our profile. It is well-known that 
this coupling is most efficient when the incident angle 
satisfies 



for Y < 1, or 



sin 9 = ± 



± 



Y 



Y + l 



■ sin a 



Y 



Y - 1 



sin a 



(26) 



(27) 



for Y > 1, where a is the angle between the external 
magnetic field and the direction of inhomogeneity and is 
equal to 90° in our casem££22£&2& Then it is easy to 
see that there is no incident angle satisfying Eq. (|2*7jl and 
that Eq. I|26|) can be rewritten as 



(28) 



v/i + w 2 /C 2/3 



where we have kept only the plus sign. 

In Fig.EJa), we plot fi, ma x and H2,max, which are the v 
values where Ai and A2 take the maximum values Ai )max 
and A 2ima x respectively in Fig. 0] versus w for ( = 500, 
v = 10~ 8 and L/A = 20 and compare them with Eq. J2HJ>- 
We find that when w is greater than about 1, the agree- 
ment between our result and Eq. I|28|l is excellent. When 
W < 1, however, our result differs greatly from Eq. 128|l . 
In Fig- Efb), the maximum mode conversion coefficients 
^4i, max and ^2, max are plotted versus w. 

In the rest of this subsection, we consider the spatial 
dependence of electric and magnetic field intensities in- 
side the inhomogeneous plasma. Using Eqs. (|14|) and 
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(|19fl . we have calculated the electric and magnetic field 
components E y i(z) and B y i{z) for incident s waves and 
E y2 (z) and B y2 (z) for incident p waves. In Fig. we 
plot the absolute values of E y i(z), B y i(z), E y2 (z) and 
B y2 {z) for = 50, £> = 10~ 3 and L/A = 20. v is chosen 
to be equal to the critical value v c = w/y^l + w 2 /£ 2 / 3 
and the value of w is 1 in Figs. E^a-d), 2 in Figs. 0Vh) 
and 3 in Figs. EJi-l). On these figures, we also indicate 
the positions of the resonant point zr given by Eq. I|25|) 
and the point z c where incident waves start to become 
evanescent when Eq. l(2*5|) is satisfied. z c turns out to be 
given by 

2Y 

z c -z Q = A^-y, (29) 

where Y = u> 2 /£ 2 / 3 . For zr < z < z c , electromagnetic 
waves are expected to be evanescent. The highly oscil- 
lating behavior for z < zr in Figs. EIi-1) comes from a 
branch of the X mode commonly known as the Z mode. 

In Fig. |S| we compare the field distributions at the 
critical value of v given by Eq. (|28|l with those for v 
slightly different from the critical value. In Figs. Ete-h), 
C = 50, v = 10~ 3 , L/A = 20, to = 3 and v = v c w 2.326. 
v is equal to v c — 0.1 in Figs. ISJa-d) and v c + 0.1 in 
Figs. Eti-l) an d all other parameters are the same as in 
Figs. ISfe-h). We find that the evanescent behavior is 
absent in off-critical cases. 

In Fig. El we show the absolute values of the z compo- 



nents of the electric field, E z \ for incident s waves and 
E z2 for incident p waves, calculated using Eq. (|A2(I . The 
parameters used are C = 50, v = 10~ 3 , L/A = 20, w = 3 
and v — v c for Figs.[?J{c-d), v = v c — 0.3 for Figs. EJa-b), 
v = v c + 0.3 for Figs. |H[e-f). In off-critical cases, the 
evanescent behavior is absent and the divergence at the 
resonance point becomes weaker. 



C. Frequency dependence 

For the study of the frequency dependence of mode 
conversion phenomena, we use the plasma density profile 
given by Eq. (|22|l . We assume /o = 1 MHz, so the local 
plasma frequency varies from 1 MHz at z = L to 5 MHz 
at z = 0. The cyclotron frequency f c is 1.5 MHz and the 
collision frequency f v is equal to 10~ 8 MHz. We show 
the mode conversion coefficients A± and A 2 for several 
incident angles and for L — 250 m in Fig. ^| and for 
L = 1500 m in Fig.^] In all cases, mode conversion is 
found to occur only in the frequency range determined 
by the upper hybrid frequency. In other words, A\ and 
A 2 are nonzero only for \/l 2 + 1.5 2 « 1.8 MHz < f < 
V5 2 + 1.5 2 w 5.22 MHz. 

In the case of uniform magnetized plasmas, the O (X) 
wave shows a resonance when the wave frequency u) is 
equal to the resonance frequency too {^x ) ■ The frequen- 
cies loo and cox are given by 



too 



\ 



LOpLO 2 COS 2 l 



\ 



L0l+L0 2 c 



(30) 



where <p is the angle between the external magnetic field 
and the wave vector^ The fact that mode conversion oc- 
curs only when lo is equal to Lu u h regardless of the incident 
angle suggests that in our stratified plasma, only the X 
wave component propagating in the direction of inhomo- 
geneity with <fi = 90° can cause mode conversion. 



V. CONCLUSION 

In this paper, we have presented a new version of the 
invariant imbedding theory for the propagation of cou- 
pled waves in inhomogeneous media and applied it to 
the mode conversion of high frequency electromagnetic 
waves into electrostatic modes in cold, magnetized and 
stratified plasmas. We have considered the cases where 



the external magnetic field is applied perpendicularly to 
the direction of inhomogeneity and the electron density 
profile is linear. We have obtained extensive and numer- 
ically exact results for the mode conversion coefficients, 
the reflectances and the wave electric and magnetic field 
profiles inside the inhomogeneous plasma. We have ex- 
plored the dependences of mode conversion phenomena 
on the size of the external magnetic field, the incident 
angle and the wave frequency in detail. Our theoretical 
method and results are expected to be quite useful in 
investigating a broad range of observations and experi- 
ments on the interaction between electromagnetic waves 
and plasmas both in space and laboratories. In forthcom- 
ing papers, we will apply our theory to the cases where 
the external magnetic field is applied in different direc- 
tions and to the propagation and mode conversion of low 
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frequency waves in more general plasmas consisting of 
various kinds of ions as well as electrons. 



of E y and B y : 
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APPENDIX A: DERIVATION OF THE 
MATRICES £ AND M 



Using the Maxwell's equations 

1 <9B 

V x E = — = ikB, 

c dt 



„ „ 1<9D 
VxB=-— 

c at 



-ike ■ E 



(Al) 



and the dielectric tensor, Eq. Ijl5(l . we express the z- 
dependent field amplitudes E x , E z , B x and B z in terms 



E x {z) 
E z {z) 



ke 3 (z^ " 
.e 2 (z) 



By'(z), 



£l(z) 



, „ sin 8 , „ 

Ey(z) - —-By(z), 



ei(» 



B x (z) = -Ey'(z), 
B z (z) = E y {z)sme, 



(A2) 



where the prime denotes a differentiation with respect to 
z. We substitute these equations and their derivatives 
into the z components of Eq. i|18|) and obtain a matrix 
equation of the form 



V>" + Aip' + B*p = 0, 



(A3) 



where 



_ ( Ey(z) 
By(z) 



B 



A = 

e 2 (z)e 3 (z) 





-■ 





53 W 



-i ) { sin 



£3(2) 




(A4) 



Comparing this with Eq. (JSJ), we get A = —£'£ 1 and / e\ — s sin 2 9 —i—sm9\ 

B = £K 2 M = k 2 £M with K = kl. From the forms of M = I gin6 , x "' sin 2 fl ) ^ A5 ) 

A and we deduce the matrices ^ Cl «i / 
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FIG. 2: Mode conversion coefficient A\ (A2) for obliquely- 
incident s (p) waves as a function of w for several values of 
v = C 1/3 sin6) and for C = 500, v = 10~ 8 and L/A = 20. 




FIG. 1: Mode conversion coefficient A\ for normally-incident 
s waves as a function of the parameter w = £ 1//3 \/Y for several 
values of C = ^A/c and for v = 10~ 8 and L/A = 20. Ai 
approaches a universal curve as £ becomes large. 




FIG. 4: Mode conversion coefficients Ai and A2 as a function 
of v for several values of w and for £ = 500, v = 10~ 8 and 
L/A = 20. 
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FIG. 5: Reflectances Rn, R22 and R12 (= R21) as a function 
of v for C = 500, v = 10" 8 , L/A = 20 and (a) w = 0.5, (b) 
w = 1, (c) w = 2. 




FIG. 6: (a) ui, ma x («2,max), which is the v value where Ai 
(A2) takes the maximum value j4i, max (^.max) in Fig.|lJ as a 
function of w for £ = 500, v = 10 -8 and L/A = 20. Our exact 
results are compared with Eq. <|280 . (b) j4i, max and ^2, max 
versus w. 
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FIG. 7: Spatial dependence of the absolute values of the elec- 
tric and magnetic field components E y i(z) and B y i(z) for 
incident s waves and E y 2(z) and B y 2(z) for incident p waves, 
for C = 50, v = 1CT 3 , L/A = 20, v = v c = w/y/l + ^ 2 /C 2/3 
and (a-d) w = 1, (e-h) w = 2, (i-1) w = 3. zr and z c are 
defined by Eqs. 1251 and 129H respectively. For zr < z < z c , 
the waves are expected to be evanescent. 
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FIG. 9: Spatial dependence of the absolute values of the z 
components of the electric field for £ = 50, v = 10 -3 , L/A — 
20, w = 3 and (a-b) v = v c — 0.3, (c-d) « = v c , (e-f) n = 
v c + 0.3. 




FIG. 10: Frequency dependence of the mode conversion coeffi- 
cients Ai and A2 for the linear profile given by Eq. I2'2ti and for 
incident angles 9 = 20°, 30°, 40°, 50° and 60°. / = 1 MHZ, 
f c = 1.5 MHZ, /„ = 10" 8 MHz and L = 250 m. Note that 
mode conversion is limited to the frequency range given by 
the upper hybrid frequency, that is \/l 2 + 1.5 2 ~ 1.8 MHz < 
/ < V5 2 + 1.5 2 « 5.22 MHz. 
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FIG. 11: Frequency dependence of the mode conversion coef- 
ficients A\ and A2 for the linear profile given by Eq. ^ and 
for incident angles 9 = 20°, 30°, 40°, 50° and 60°. /o = 1 
MHZ, f c = 1.5 MHZ, f v = 10" 8 MHz and L = 1500 m. 



